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Confexf. Bouncing collisions of dust aggregates within the protoplanetary may have a significant impact on the growth process of 
planetesimals. Yet, the conditions that result in bouncing are not very well understood. Existing simulations studying the bouncing 
behavior used aggregates with an artificial, very regular internal structure. 

Aims. Here, we study the bouncing behavior of sub-mm dust aggregates that are constructed applying different sample preparation 
methods. We analyze how the internal structure of the aggregate alters the collisional outcome and determine the influence of 
aggregate size, porosity, collision velocity, and impact parameter. 

Methods. We use molecular dynamics simulations where the individual aggregates are treated as spheres that are made up of several 
hundred thousand individual monomers. The simulations are run on GPUs. 

Results. Statistical bulk properties and thus bouncing behavior of sub-mm dust aggregates depend heavily on the preparation method. 
In particular, there is no unique relation between the average volume filling factor and the coordination number of the aggregate. 
Realistic aggregates bounce only if their volume filling factor exceeds 0.5 and collision velocities are below 0.1 ms -1 . 
Conclusions. For dust particles in the protoplanetary nebula we suggest that the bouncing barrier may not be such a strong handicap 
in the growth phase of dust agglomerates, at least in the size range of « 100/im. 
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1. Introduction 

For the planet formation process, the growth from micron sized 
dust grains to kilometer sized objects is a key ingredient of the 
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J core accretion scenario originally proposed by |Pollack et al. 
( 1996). Yet, the question remains unanswered how this process 
is accomplished in the face of various impediments. First of all, 
fast inward drift limits the time available to for m planetesimals 
by successive collisions to less than 10 4 years ( Weidenschilling 
1977}, The growth process itself heavily depends on two ingre- 



dients: 

1 . The dynamical properties of the disk that determine the col- 
lision rate as well as the parameters of a collision such as the 
impact velocity. 

2. The mechanical behavior of the colliding aggregates that de- 
termines the outcome of the collision. 

Since the information accessible through direct observations is 
limited the first aspect is addressed mainly by theoretical work 
and computer simulations (Flaig e t al.|2012| . For an overview on 
the properties of protoplanetary disks we refer to the following 
reviews by |Dullemond et aL ( |2007| l and | Armitage| ( |20 IT) . 

To investigate the collision behavior of dust/icy aggregates 
various methods are employed. Depending on the size of the 
aggregates and the desired collision velocity, laboratory exper- 
iments are possible. As of today, laboratory experiments pro- 
vide data of collisions ranging from millimeter- to decimeter- 
sized aggregates composed of different materials (mainly Sili- 
cate/Quartz and Ice). A comprehensive summary of laboratory 
experiments is given by Blum & Wurm (2008 ). Computer simu- 
lations provide a second method to study the collisional behavior 



of dust or ice aggregates. Collisions of very small, micron sized 
aggregates have been simulated using a molecular dynamics ap- 
proach (e.g. Dominik & Tielens|1997||Paszun & Dominik[ 2009 ; 
Wada et al. 2007 2009). For macroscopic aggregates different 



methods such as smoothed particle hydrodynamics (SPH) are 
employed (e.g. Sch afer et al.|2007]|Geretshauser et al.|2010 1 



More recent experiments showed that collisions of mm to 
cm-sized aggregates often result in bouncing (e.g. 
~E009 " 



et all 



Weidling 

|HeiBelmann et al.|[20T0l |Weidling et al.||2012[ 



Jankowski et al.|2012| ). Extrapolating the results obtained from 



the various experiments Guttler et al.| ( |2010| ) devised a model 
describing the outcome of collision with respect to the collision 
velocity, and the mass and porosity of the colliding aggregates. 
Employing this model to simulate the evolution of a swarm of 
dust aggregates in a protoplanetary disk the so called "bouncing 
barrier" emerged (Zsom et al. 20T0|). As the aggregates grow 
larger their relative velocities increase. Due to the growing ki- 
netic impact energy aggregates get increasingly compacted dur- 
ing successive collisions. When the aggregates get too compact 
their collisions do not result in sticking anymore. Instead, they 
bounce off each other and the growth process is stopped. This 
occurs in the size regime of centimeters. 

A possible way to overcome the bouncing barrier has been 
recently suggested by Windmark et al. (2012a). Under the as- 
sumption of a few bigger aggregates that act as initial seeds it 
is possible to grow larger 100 m sized bodies on the timescale 
of 1 Myr. A possible origin of those seeds has been proposed 
by |Windmark et aL] ( |2012b| ). Taking into account a Maxwellian 
velocity distribution they found that low velocity collisions can 
allow a few aggregates to grow considerably larger than the av- 
erage of the simulated population. 
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Despite its significant influence on the growth process, 
bouncing still lacks theoretical understanding of its prerequisites 
on a microscopical scale. So far, molecular dynamics (MD) sim- 
ulations result in bouncing only for rather compact aggregates 
(|Wada et al.|[20TT) |Schrapler et"aLl|20"T2"l >. According to |Wada 



Table 1. Material Parameters of the individual monomers used in the 
simulations. 



et al. ( 201 1) an average coordination number of 6 is required for 
aggregates to bounce off each other. However, in laboratory ex- 
periments bouncing frequently observed in collisions of aggre- 
gates with lower filling factors / coordination numbers for which 
MD simulations clearly predict sticking. It has been speculated 
that this discrepancy could result from a size effect or a possible 
compaction of the outer shell during the handling process of the 
aggregates used in the laboratory experiments. However, the lat- 
ter hypothesis has been refuted by recent experiments (Kothe et 
al. 2012). 

The aim of this work is to study the influence of the internal 
structure on the bouncing behavior of sub-mm dust aggregates in 
greater detail. Using the enormous computing power provided 
by GPUs we can simulate aggregates consisting of several hun- 
dreds of thousands of monomers and thus simulate aggregates in 
a size range from several microns up to « 0.1 mm in diameter. 

2. Interaction model 

To simulate the behavior of of dust aggregates we use a soft 
sphere discrete element method (SSDEM) approach. The dust 
aggregates are composed of hundreds of thousands of micron 
sized spherical grains (monomers). Our interaction model is 
based on the work of [Dominik & Tielens (1997]) who com 
bined earlier theoretical work by |Johnson et al.| ( 1971 
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mini k & Tielens| ( |1995| |1996) into a detailed micro-mechanical 
model describing the interaction between two monomers. These 
monomers may establish adhesive contacts when touching each 
other and kinetic energy is dissipated upon deformation of these 
contacts. A few years later, |Wada e t al. (2007 ) presented a differ- 
ent approach when deriving nearly the same forces and torques 
from corresponding potentials. 

When trying to reproduce the results of laboratory experi- 
ments performed b y|Giittler et al.|([2009|l on the compression of 
porous dust cakes Seizinger et al. |(|2012[) observed t hat the be- 
havior predicted by the model of |Dominik & Tielens ( 1997} was 
too soft. Since the samples used by Guttler et al.| ( 20090 had 
been composed of micron-sized, spherical, monodisperse sili- 
cate grains their results constituted a perfect possibility to cali- 
brate the model. Introducing two free parameters m r and m s that 
modify the strength of the rolling and sliding interaction between 
two monomers Seizinger et al. ( 2012| > were able to obtain excel- 
lent agreement between laboratory results and computer simula- 
tions. 

In this work we use the modified interaction model presented 

by 



Seizinger et al. 



parameters are 



2012| ) with m r — 8 and m & — 2.5. The mate- 



isted in Tab.[T] 



3. Sample generation 

In this work we examine the conditions under which bounc- 
ing occurs. Apart from the external parameters describing the 
physics of collisions such as the impact parameter or velocity 
we study the influence of the internal structure of the aggregate. 
Examples of such aggregates that have been generated by differ- 
ent methods are shown in Fig.[T] 

To study the influence of the aggregate size we use aggre- 
gates with diameters in the range of 30 to lOOyi/m. Unfortu- 
nately, simulations with larger aggregates are infeasible due to 



Physical property 


Silicate 


Particle Radius r (in /mi) 


0.6 


Density p (in g cm -3 ) 


2.65 


Surface Energy y (in mJ irT 2 ) 


20 


Young's Modulus E (in GPa) 


54 


Poisson Number v 


0.17 


Critical Rolling Length £ cl ; t (in nm) 


2 



the required computational cost, at least for a wider range of pa- 
rameters. 

Typically, the volume filling factor (f> and the average coordi- 
nation number « c are used to classify aggregates. In general, the 
filling factor is given by 



NV P 
V A 



(1) 



where N denotes the number of monomers, V p is the volume of a 
monomer, and Va is the total volume occupied by the aggregate. 
As we use spherical aggregates Va can be calculated easily from 
the outer radius of the aggregate. For irregular shaped aggre- 
gates there are different possibilities to define Va- For example 
|Kozasa et al.| ( |1992| l determine the size of a porous aggregate 
based on its radius of gyration whereas Ossenkopf (1993) use 
the geometric cross section. It is important to keep this ambigu- 
ity in mind when comparing the porosities of fluffy aggregates 
to other works. 

In molecular dynamics simulations, the coordination number 
of a monomer denotes the number of the neighbors it interacts 
with. Thus, the average coordination number « c is obtained by 
averaging the number of contacts of each particle. 

In this work we use three different types of aggregates: 
Hexagonal lattice aggregates featuring a regular grid struc- 
ture, aggregates produced by slowly compacting a porous dust 
cake, and aggregates generated by successively adding single 
monomers from randomly chosen directions. These choices 
have been motivated by the fact that hexagonal lattice aggregates 
are easy to build and allow for comparison with earlier work by 
Wada et aX](|201 1), whereas the static compaction resembles the 



generation of samples used in laboratory results. The aggregates 
of the third type are generated algorithmically but their structure 
remains comparable to the static compaction type (see Sect. 3.4 1. 

Because the aggregates within the protoplanetary nebula 
grow through successive collisions one might expect that their 
internal structure lies somewhere in between the static compres- 
sion and the ballistic aggregation cases. 

3.1. Hexagonal lattice (CPE) 

Hexagonal-lattice type aggregates (also referred to as hexago- 
nal close packing with extraction (CPE)) may be generated very 
easily. First, a hexagonal close packing aggregate is generated 
which features a volume filling factor <p ~ 0.74 and a coordina- 
tion number n c m 12 (due to surface effects n c equals 12 only for 
aggregates of infinite size). In the second step a suitable num- 
ber of randomly selected monomers will be removed to achieve 
the desired volume filling factor. As a result a small number of 
monomers on the surface may be become disconnected from the 
main aggregate and will be removed as well. 
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Fig. 1. Examples of the different types of aggregates used in this work: (a) Hexagonal lattice <f> = 0.59, n c = 9.93, (b) Ballistic aggregation with 
migration (f> = 0.40, n c = 3.98, and (c) Static compaction = 0.49, n c = 3.50. All depicted spheres have a diameter of 60 yum. 



Wada et al. ( 201 l| l have already studied the bouncing behav- 
ior of this type of aggregates and found that bouncing will oc- 



cur if the average coordination number is greater than 6. Schra- 



pler et al. ( 2012| l analyzed the relation between the coefficient 



of restitution and the collision velocity in experiments and sim- 
ulations using CPE-aggregates. Their results agreed well with a 
theoretical model by Thornton & Ning ( 1998 ). 



3.2. Ballistic aggregation with migration (BAM) 



The second type of aggregates was originally suggested by Shen| 
[eTaLl ( [2008] l and also studied in the work of | Wada et aT] l |2011) . 
To generate a larger aggregate single monomers are succes- 
sively shot in from random directions onto the existing aggre- 
gate. When the monomer hits the aggregate it will either remain 
at the position where the first contact has been established or 
migrate to a position close by where it establishes contacts with 
two or three monomers. Compared to Shen et al. (2008 1, we use 
three different methods to select the final position of the migrat- 
ing particle: 

1. Select the position closest to the spot, where the monomer 
impacts on the aggregate (referred to as "shortest migra- 
tion"). 

2. Select the position randomly from all available possibilities 
(referred to as "random migration"). 

3. Select the position which is closest to the center of mass (re- 
ferred to as "center migration"). 

For a given coordination number the resulting aggregates show 
a different filling factor depending on which selection mecha- 
nism is employed (see Fig. [2}. The first method leads to rather 
porous aggregates since the monomers typically migrate to posi- 
tions further outward compared to the case of random migration. 
Likewise, the resulting aggregates will become even more com- 
pact if monomers migrate to the most inward position available. 

Note that we do not claim that random or center migration 
are realistic growth processes that accurately describe the growth 
of dust aggregates in protoplanetary discs. Yet, they constitute 
a computationally very cheap approach to generate larger aggre- 
gates that do not suffer from the artificial lattice structure like 
the CPE aggregates described above. Compared to the "static 
compaction"-aggregates they additionally offer the advantage 
that they are perfectly relaxated. Since all monomers are in equi- 
librium distance from each other, there are no attractive or repul- 
sive forces that could lead to a breakup of the aggregate. 
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Fig. 2. The relation between the volume filling factor <p and the av- 
erage coordination number n c . All aggregates are spherical and have a 
diameter of 60 jim. 



3.3. Static compaction (SC) 

The last type of aggregates used for our the studies is the 
most computationally expensive. To generate a spherical ag- 
gregate of a certain diameter we start with a sufficiently large, 
cuboid shaped dust cake generated by random ballistic depo- 
sition (RBD). Since RBD-aggregates feature an initial volume 
filling factor of <p — 0.15 we first have to compact the aggregate 
until we reach the desired filling factor. For this purpose, the 
aggregate is put into a box of walls that may move towards each 
other. According to Seizinger et al. (2012 1 this compaction must 
be very slow to avoid inhomogeneities. 

Even if the cake is compacted homogeneously for filling fac- 
tors above m 0.45 it will get elastically charged and expand when 
the compacting walls are removed. Therefore the dust cakes 
needs to be relaxated before removing the walls of the com- 
paction box. For this purpose we disturb the aggregate by mod- 
ifying the position of the monomers randomly by a very small 
amount. We keep the disturbed aggregate in a box of fixed size 
and wait until the kinetic energy induced by the disturbance is 
damped away by the inelastic monomer interaction. To get rid 
of kinetic energy below the threshold where the inelastic regime 
is entered we additionally enforce a viscous damping mecha- 
nism. For this purpose, the velocities and angular velocities of 
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the monomers are multiplied by a factor of 1 - k in each time 
step, where k denotes a damping coefficient. In this work we use 
k = 0.0001. 

It turns out that a maximum disturbance of a factor of 0.001 
of the radius of a monomer is sufficient to stabilize the aggregate 
without altering its internal structure. Higher values may alter 
the coordination number significantly which could potentially 
change the collisional behavior of the aggregates and is therefore 
unwanted. For fillings factors above <p = 0.58 - 0.59 this method 
does not work anymore. Here the compaction is too close to the 
random closest packing (cp as 0.63) and no stable configuration 
can be reached without rearranging the monomers significantly. 

After the aggregate has been relaxated the compaction box 
is removed and a spherical aggregate will be cut out of the com- 
pacted cake. As this procedure is computationally very expen- 
sive it takes several days to generate larger (50 fim in diameter 
and above) aggregates of this type. 





3.4. Comparison 

Comparing the relation between and n c of the different types 
of aggregates described above shows very interesting features: 
As we can see in Fig. [2] the relation depends considerably on 
the preparation method. The different BAM generation methods 
have been described in Sect. 13. 21 

From the work presented in this section two important con- 
clusions can be drawn: 

1 . The coordination number n c is not sufficient to describe the 
properties of an aggregate. Especially, there is no unam- 
biguous relation (p(n c ) between filling factor and coordina- 
tion number. 

2. Hexagonal lattice (CPE) aggregates have a very distinct re- 
lation <p{n c ) compared to the other two methods that produce 
aggregates with less artificial structures. 

In laboratory experiments with aggregates composed of micron 
sized dust grains, it is typically only possible to determine the 
filling factor but not the coordination number. Thus, one has to 
be very careful when comparing results from numerical simula- 
tions of CPE aggregates with laboratory experiments. 

4. Results 

In this section we present our results from various simulations 
in which we study the influence of the collisions velocity, im- 
pact parameter, and aggregate size on the bouncing behavior. All 
simulations have been performed on NVIDIA GPUs (GTX460, 
GTX570, Tesla C2070). Depending on the aggregate size and 
filling factor each simulation took between less than an hour and 
half a day. 

4. 1 . Growth Factor 

In the following bouncing maps the "growth factor" y that is in- 
spired by the four-population model suggested by |Geretshauser| 



et al. (201 1 1 is depicted. It is defined by 



r : 



"^largest 
'«tot 



(2) 



where mi argest is the mass of the largest fragment and m tot the total 
mass of the colliding aggregates. For perfect sticking we obtain 
7=1, for total destruction y — > 0. In collisions of equal sized 
aggregates, a value y = 0.5 indicates bouncing. However, during 




Fig. 3. Schematic view of the possible alignment of two CPE aggre- 
gates. In the upper case the aggregates are oriented with respect to their 
lattice structure. In the lower case the orientation is arbitrary. 



the transition from perfect sticking to fragmentation y may also 
become 0.5. To distinguish between the two cases we consider 
the mass ratio 72 of the second largest fragment. In the bouncing 
case it is 0.5 as well whereas in the fragmentation case the mass 
of the second largest fragment is much lower than 0.5 of the total 
mass as there are a lot of other smaller fragments. 

Thus, in the bouncing maps presented in this work the green 
areas indicate sticking, the upper left yellow areas bouncing, and 
the color gradient from green to yellow to red on the right marks 
the transition from sticking to fragmentation. 

Note that y - yi - 0.5 only applies in the case of "per- 
fect bouncing". In our simulations we often observe the loss 
of a few monomers (typically less than 100) which is negligible 
compared to the total number of monomers of 5 • 10 4 to 5 • 10 5 . 
Thus, we also count collisions as bouncing events if y and yi are 
slightly smaller than 0.5. 



4.2. Hexagonal lattice 

The outcome of head-on collisions of CPE aggregates has al- 
ready been studied by| Wada et al. ( 201 1 ) who observed bouncing 



if the coordination number was greater than 6. However, their 
aggregates were much smaller (w 10 4 monomers). As hexagonal 
lattice aggregates feature a regular lattice structure their orienta- 
tion is likely to influence the collision behavior. Thus, we first 
examine the effect the orientation by comparing the case where 
the aggregates are aligned to their lattice structure (see upper 
part of Pig.[3j to a random orientation (lower part of Fig. [3}). 
As we can see in Fig.|4 the orientation of the aggregates is im- 
portant especially for the transition from sticking to bouncing 
with increasing filling factor. Looking at the left panel of Fig. [4] 
and comparing the filling factor with the coordination number in 
Fig.|2] we can reproduce the n c > 6-criterion proposed by |Wada| 
|et al. ( |201 \\ for the aligned case. On the other hand, the bounc- 
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Fig. 4. Growth factor, y (eq. [2|, of the collision of two CPE aggregates with a diameter of 60 fim. Sticking occurs in the green colored area 
whereas the yellow area in the upper left indicates bouncing. Left: Orientation aligned to the lattice structure of the aggregates. Right: Non 
aligned orientation. 
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Fig. 5. Growth factor of the collision of two CPE aggregates of different size averaged over three different orientations. Left: Aggregates of a 
diameter of 30 /im. Right: Aggregates of a diameter of 60 //m. 



ing maps differs significantly for a non aligned orientation (see 
right panel of Fig. [4]). 

In order to mitigate the effect of the orientation we averaged 
over three different orientations to investigate the size depen- 
dency of our results. Each map has been generated using 12 
different filling factors and 28 velocities. Thus, 3 • 336 = 1008 
simulations had to be performed in total. Concerning bouncing 
we could not observe a clear difference between aggregates with 
a diameter of 30 and 60 yum (see Fig. [5]). 

However, the velocity v s ^ t at which the transition from stick- 
ing to fragmentation occurs changes significantly. For the small 
aggregates we get v s ^f * 4ms~' (left panel of Fig. [5]). For 
the bigger aggregates we observe that v s -tf depends on the fill- 
ing factor. For (f> < 0.43 we get v s _,f ~ 10ms _I whereas 
v s _>f ~ 5ms~' for <p > 0.43 (right panel of Fig. [5}. This can 
be explained by the reduced capability of compact aggregates to 
dissipate kinetic energy by restructuring. Taking into account 
Fig.|2]we see that the transition occurs when the average coordi- 
nation number « c exceeds a value of 6. A monomer with six or 
more contacts is fixated rather tightly and and thus the aggregate 
cannot change its internal structure as easily anymore. 



In summary it can be said that for hexagonal lattice aggre- 
gates we regularly observe bouncing collisions for filling factors 
above 0.5 and collision velocities up to roughly 0.2ms -1 . 

4.3. Ballistic aggregation with migration 

Wad a et aT] ( |201 1[ ) found that bouncing may occur if n c > 6 
independent of the type of aggregate they used. In Fig. [6] we 
show the outcome of collisions between two roughly 75 fim 
sized BAM aggregates generated by using the shortest migration 
method described in Sect. |3.2| The corresponding filling factor 
is between 0.36 and 0.39. However, we did observe only two 
bouncing collisions. Since n c = 6 is the maximum value that 
can be achieved by two times migration we could not investigate 
what happens at higher coordination numbers. 

Repeating the setup described above for the center migra- 
tion case we get similar results as for the shortest migration case 
shown in Fig. [6] We observe hardly any bouncing events even 
for aggregates with n c = 6 (which corresponds to a filling fac- 
tor of 0.49 - 0.5). This indicates that the bouncing behavior 
of BAM aggregates depends more on the filling factor than the 
coordination number. 
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Collision Speed v [m/s] 

Fig. 6. Growth factor of the collision of two BAM aggregates gen- 
erated by the shortest migration method. For three different aggregates 
of a diameter of roughly 75 /jm only two collisions at 1 cms -1 lead to 
bouncing. 



To achieve higher filling factors we switched to the center 
migration method (see Sect. |3.2[ ). The corresponding bouncing 
maps are shown in Fig. [7] It is striking that the bouncing regime 
is much smaller compared to the CPE aggregates. As before, for 
larger aggregates the transition from sticking to fragmentation 
occurs at higher velocities. 

4.4. Static compaction 

The bouncing behavior of the aggregates generated by static 
compaction is quite similar to the BAM aggregates (see Fig. [8}. 
Again, the bouncing regime is considerably smaller compared 
to the case of hexagonal lattice aggregates and bouncing is ob- 
served only in some cases for high filling factors above 0.5 and 
collision velocities below 0.1 ms _1 . 

As in the case of the other aggregate types the transition 
velocity v s ^f from sticking to fragmentation increases with in- 
creasing aggregate size. For the the small aggregates with 
d = 30//m we observe v s ^f m 4ms~ I whereas for d = 60 //m 
the transition velocity goes up to v s _,f ~ 12 ms _1 . 

4.5. Size dependency 

To further examine the influence of the aggregate size we per- 
formed collisions of 100//m-sized CPE and BAM aggregates 
(using the center migration method as it yields the BAM aggre- 
gates with the highest filling factors). For the CPE aggregates 
we observe slightly more bouncing for filling factors between 
0.4 to 0.5 (see left panel of Fig.|9]l. However, for BAM aggre- 
gates there is no noticeable difference compared to the 60 //m 
aggregates (see right panel of Fig.[9]l. 

Depending on the filling factor the lOO/vm aggregates con- 
sist of up to 350, 000 monomers. In order to analyze the size 
dependency it would be desirable to simulate collisions of even 
larger aggregates. Unfortunately, this is rendered impossible by 
the lack of available computing power. Simulating a single colli- 
sion of two 100 fim sized aggregates took 10 to 20 hours (due to 
the different filling factors) on a GPU. Doubling the size would 
require computing times on the order of weeks for a single col- 
lision. For each orientation shown in Fig. [9] 48 collisions have 
been performed. 



4.6. Impact parameter 

As a last step we examine the influence of the impact parameter 
b = 0.5. For the collisions, we used the same orientation as for 
the results shown in the right panel of Fig.|4] Contrary to our 
expectations we do not observe a significant influence of the im- 
pact parameter on the bouncing behavior of CPE aggregates (see 
Fig. [TO}. However, fragmentation sets in at considerably lower 
velocities of v s ^f * 3ms _1 . In a head on collision the entire 
aggregate may dissipate the kinetic impact energy by internal re- 
structuring and thus help to avoid fragmentation. This does not 
apply to offset collisions where it is easier to tear away the outer 
layers without major restructuring of the core of the aggregates. 

However, comparing the growth map of the collisions be- 
tween aggregates with a diameter of 30 and 60 pun we do not 
observe any significant increase of velocity v s ^f where the tran- 
sition from sticking to fragmentation occurs (see left and right 
panel of Fig.fTO}. As already pointed out by |Wada et al.| (2009) 
the increase of v s _,f for larger aggregates is limited to the case of 
head-on collisions. 

In contrast, for BAM aggregates we do observe bouncing 
in a larger regime compared to case of head-on collisions (see 
Fig.[TT|l. Again, the same samples and orientation as for the 
head-on collisions (right panel of Fig.[7]i have been used. 



5. Requirements for bouncing 

Compared to the aggregates generated by BAM or static com- 
paction the bouncing regime of CPE aggregates is significantly 
larger. A likely explanation for this discrepancy is given by the 
different structure of the aggregates. To gain a deeper insight 
in the physical processes leading to sticking or bouncing it is 
worthwhile to have a closer look at a single collision. 

Two aggregates may bounce off each other only if there is 
enough elastic energy left to break the contact area. Thus, a 
significant amount of the kinetic impact energy must be stored 
temporarily without being dissipated. If the colliding aggregates 
penetrate each other too deeply the impact energy is dissipated 
upon internal restructuring in the area where the contact is es- 
tablished. In this case inelastic sliding and rolling constitute the 
main dissipation channels ( |Wada et al.|201 \\ . Thus, the ratio of 
elastic to dissipated energy of colliding aggregates is the key pa- 
rameter that determines whether sticking or bouncing will occur. 

Being able to track the evolution of the different types of 
energies over time is the key advantage of the model presented 
by |Wada et"aE] ( |2007) l. To address the different behavior of BAM 
and CPE aggregates we compare a bouncing collision of two 
CPE aggregates with a sticking collision of BAM aggregates. 
Both aggregates are 60 //m in diameter and have a filling factor 
<p ~ 0.59. The time evolution of different types of energies and 
potentials for such collisions is shown in Fig. [12] 

As expected, in the sticking case most kinetic energy is dis- 
sipated by inelastic sliding and rolling (right panel of Fig.[T2|. 
Only a small percentage of the impact energy is stored in the 
elastic regime of the normal f/ norma i and sliding potential E/ S iide 
(since the elastic energy stored in the rolling and twisting poten- 
tials is negligible they are not shown in Fig. 12 1. 

We observe an entirely different situation in the bouncing 
case: As shown in the left panel of Fig. 12 only about one third 
of the impact energy is dissipated, whilst roughly half of the ki- 
netic energy is temporarily converted into potential energy. This 
coincides well with our predictions above. The potential energy 
that is stored mainly in the normal and sliding interaction is con- 
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Fig. 8. Growth factor of the collision of two static compaction aggregates of different size. Left: Diameter of 30/jm. Right: Diameter of 60/um. 



verted back into kinetic energy and allows that the colliding ag- 
gregates to separate again. 

We can conclude that due to their lattice structure CPE aggre- 
gates can convert significantly more impact energy into elastic 
energy than BAM aggregates. In a compact CPE aggregate the 
monomers are located in densely packed layers. When the outer 
monomer of such a layer hits the other aggregate it is pushed in- 
wards and will compress the layer. This way, kinetic energy is 
converted into potential energy without the occurrence of inelas- 
tic restructuring. This mechanism works well in the presence 
of a regular grid structure as it is the case for CPE aggregates. 
However, the monomers of BAM aggregates are not arranged in 
any regular pattern. Thus, they are not likely to bounce unless 
they are very compact in which case energy dissipation by inter- 
nal restructuring is hindered because the monomers are locked 
in their position. 

It also offers an explanation for the lower impact velocity of 
BAM aggregates at which we observe the transition from bounc- 
ing to sticking. At impact velocities above 0.1 ms~' the kinetic 
impact energy is sufficient to restructure the monomers in the 
contact area that had been locked at lower impact velocities. 
Since the lattice structure of CPE aggregates offers higher resis- 
tivity against restructuring their transition velocity from bounc- 
ing to sticking is roughly 0.3 ms _I . 



6. Conclusions 

From the analysis of the statistical properties of the different 
samples presented in Sect. 3] we clearly see that the preparation 
method plays a crucial ro e when studying the collisional be- 
havior of microscopic dust aggregates. As the relation between 
the volume filling factor and the average coordination number 
strongly depends on the preparation method one must be care- 
ful when comparing results obtained from computer simulations 
with laboratory work. In most laboratory experiments, only the 
volume filling factor can be measured, while the likewise im- 
portant coordination number remains unknown. Since the more 
compact aggregates used in laboratory experiments are typically 
produced by mechanical compression, we expect that their mi- 
croscopic structure resembles the static compaction aggregates 
much more closely than the hexagonal lattice type aggregates. 

For computer simulations, generating large, compact aggre- 
gates by static compaction is infeasible due to the additional 
computational effort. We suggest to use BAM aggregates as an 
alternative. Their statistical properties are close to the aggre- 
gates produced by static compaction, yet they can be generated 
directly. Additionally, one does not run into problems caused by 
elastic charging as the generation procedure ensures that BAM 
aggregates are perfectly relaxed. At least in the various collisions 
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simulations performed in this work BAM and static compaction 
aggregates show very similar behavior. 

Based on the outcome of the simulations presented in this 
work we can conclude that bouncing collisions of dust aggre- 
gates in the size regime below 0.1 mm are rare. Unless the ag- 
gregates feature a regular lattice structure, which is not likely to 
be the case for the aggregates in a protoplanetary disk, bouncing 
requires filling factors greater than 0.5 and collision velocities 
below 0.1 ms~'. Even if these prerequisites are met bouncing 
does not occur very frequently. Additionally, laboratory experi- 
ments on dust growth show that the maximum filling factor that 
is achieved during the growth process is much lower than <p = 0.5 
(e.g. |Teiseret al.|2011) . Therefore, the influence of bouncing on 
the growth process is limited in the sub-mm size regime. 

Hence, from a microscopic view it remains unclear how cm- 
sized aggregates with filling factors considerably below 0.5 are 
able to bounce off each other. The idea of a compacted outer 
layer (also referred to as a hard shell) has been put forward as 
a possible explanation. Using SPH-simulations, Geretshauser et 
al. (in preparation) observed that such a hard shell can indeed 
lead to bouncing collisions between aggregates with a porous 
core. Wada et al. ( 2011| > obtained similar results when perform- 



gates featuring a hard sphere. Langkowski et al. (2008) found 
that molding an aggregate significantly alters the outcome of a 
collision experiment. However, |Kothe et al.| ( |in prep. ) analyzed 
aggregates used in their collision experiments with X-ray com- 
puter tomography imaging and could not find any compacted 
outer layers. 

In Sect. 14. 61 we have shown that offset collisions result in 
bouncing somewhat more often than head-on collisions. De- 
pending on the experimental setup head-on collisions will be 
rare, and in a free collision with many particles setup head-on 
collisions will be rare as well (e.g. Weidling et al. 2012| Beitz 
et al. 2012|l. Thus, the impact parameter helps to resolve some of 



ing molecular dynamics simulations of collisions of CPE aggre- 



the discrepancies between numerical simulations and laboratory 
experiments. 

Taking into account the different aggregate types we can only 
partially confirm the n c > 6 criterion for bouncing proposed by 
Wad a et aT] ( |201 It agrees well with our results from collisions 
of CPE aggregates. However, for BAM aggregates generated by 
the shortest migration method we observe very little bouncing at 
n c = 6 (see Sect. |4.3| and Fig.[6]i. Likewise, for static compaction 
aggregates or BAM aggregates generated by the random or cen- 
ter migration method there are a few bouncing events where n c 
is considerably lower than 6. A volume filling factor of » 0.5 
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Fig. 12. Time evolution of the kinetic, potential and dissipated energy during collisions of two aggregates with a diameter of 60 yum and filling 
factor <p « 0.59. Left: A bouncing collision between two CPE aggregates. A sufficient amount of the impact energy is temporarily stored in the 
potential (/ norma i of the normal interaction. Right: During a sticking collision of two BAM aggregate most kinetic impact energy is dissipated. 



appears to constitute a lower limit for bouncing. At the present 
time we do not have an explanation what determines the exact 
value of the critical filling factor for the onset of bouncing. We 
have shown that it depends on the amount of energy that can 
be stored in the normal potential. The energy deposition is a 
continuous process, and it is to be expected that below a certain 
density sticking ensues. Numerically, we determined this value 
to be <p « 0.5 in the sub-mm size regime. 

Our simulations give insight into the fragmentation thresh- 
old as well. For small aggregates (30 fim) the fragmentation ve- 
locity is around 4 m/s. Upon increasing the projectile size the 
fragmentation threshold rises up to about 10 m/s for the largest 
particle sizes we considered (60 /-im), and this is independent 
of the sample generation method. This value is in very good 
agreement with the findings of SPH simulation for much larger 
objects (Geretshau ser et aT]|2011| l. The shift to a larger frag- 
mentation velocity is caused by the fact that larger particles can 
dissipate more energy than smaller particles. Upon increasing 
the filling factor, the fragmentation threshold decreases because 
the aggregates become much stiffer and cannot be deformed so 



easily. However, as shown in Sect. 4.6 this effect only applies to 
the case of head-on collisions. 



With respect to the growth of small dust agglomerates in the 
protoplanetary nebula our results indicate that for more realis- 
tic aggregates (BAM-type) bouncing only occurs for very small 
collision velocities (< 0.1 m/s) and large filling factors > 0.5. 
Thus, the bouncing barrier may not be such a strong handicap in 
the growth phase of dust agglomerates, at least in the size range 
of as 100 /vm. For larger, m-sized particles SPH results indicate 
bouncing up to 1 m/s. 

Acknowledgements. A. Seizinger acknowledges the support through the German 
Research Foundation (DFG) grant KL 650/16. The authors acknowledge support 
through DFG grant KL 650/7 within the collaborative research group FOR 759 
The formation of planets. We also thank the referee, Koji Wada, for his helpful 
comments that improved the quality of this paper. A part of the simulations were 
performed on the bwGRiD cluster, which is funded by the Ministry for Education 
and Research of Germany and the Ministry for Science, Research and Arts of the 
state Baden-Wurttemberg. 



References 

Armitage, P. J. 201 1, ARA&A, 49, 195 

Beitz, E., Guttler, C, Weidling, R., & Blum, J. 2012, Icarus, 218, 701 
Blum, J. & Wurm, G. 2008, ARA&A, 46, 21 

Dominik, C. & Tielens, A. G. G. M. 1995, Philosophical Magazine, Part A, 72, 
783 



Article number, page 9 of 



10 



Dominik, C. & Helens, A. G. G. M. 1996, Philosophical Magazine, Part A, 73, 
1279 

Dominik, C. & Tielens, A. G. G. M. 1997, ApJ, 480, 647 

Dullemond, C. P., Hollenbach, D., Kamp, I., & D'Alessio, P. 2007, Protostars 
and Planets V, 555 

Flaig, M., Ruoff, P., Kley, W., & Kissmann, R. 2012, MNRAS, 420, 2419 
Geretshauser, R. J., Meru, R, Speith, R., & Kley, W. 201 1, A&A, 531, A166 
Geretshauser, R. J., Speith, R., Guttler, C, Krause, M., & Blum, J. 2010, A&A, 
513, A58 

Guttler, C, Blum, J., Zsom, A., Ormel, C. W., & Dullemond, C. P. 2010, A&A, 
513, A56 

Guttler, C, Krause, M., Geretshauser, R. J., Speith, R., & Blum, J. 2009, ApJ, 
701, 130 

HeiBelmann, D., Blum, J., Fraser, H. J., & Wolling, K. 2010, Icarus, 206, 424 

Jankowski, T., Wurm, G, Kelling, T., et al. 2012, A&A, 542, A80 

Johnson, K. L., Kendall, K, & Roberts, A. D. 1971, Royal Society of London 

Proceedings Series A, 324, 301 
Kothe, S., Blum, J., Weidling, R., & Guttler, C in prep. , submitted to Icarus 
Kozasa, T., Blum, J., & Mukai, T. 1992, A&A, 263, 423 
Langkowski, D., Teiser, J., & Blum, J. 2008, ApJ, 675, 764 
Ossenkopf, V. 1993, A&A, 280, 617 
Paszun, D. & Dominik, C 2009, A&A, 507, 1023 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62 
Schafer, C, Speith, R., & Kley, W. 2007, A&A, 470, 733 
Schrapler, R., Blum, J., Seizinger, A., & Kley, W. 2012, ApJ, 758, 35 
Seizinger, A., Speith, R., & Kley, W. 2012, A&A, 541, A59 
Shen, Y., Draine, B. T., & Johnson, E. T. 2008, ApJ, 689, 260 
Teiser, J., Engelhardt, I., & Wurm, G. 201 1, ApJ, 742, 5 
Thornton, C. & Ning, Z. 1998, Powder Technology, 99, 154 
Wada, K, Tanaka, H., Suyama, T, Kimura, H., & Yamamoto, T. 2007, ApJ, 661, 
320 

Wada, K, Tanaka, H., Suyama, T, Kimura, H., & Yamamoto, T. 2009, ApJ, 702, 
1490 

Wada, K, Tanaka, H., Suyama, T, Kimura, H., & Yamamoto, T. 201 1, ApJ, 737, 
36 

Weidenschilling, S. J. 1977, MNRAS, 180, 57 
Weidling, R., Guttler, C, & Blum, J. 2012, Icarus, 218, 688 
Weidling, R., Guttler, C, Blum, J., & Brauer, F. 2009, ApJ, 696, 2036 
Windmark, F, Birnstiel, T, Guttler, C, et al. 2012a, A&A, 540, A73 
Windmark, F, Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2012b, A&A, 
544, L16 

Zsom, A., Ormel, C. W., Guttler, C, Blum, J., & Dullemond, C. P. 2010, A&A, 
513, A57 



Article number, page 10 of 



10 



